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O ■ Abstract 

o. 

f-H , Bose-Einstein correlations are being exploited to obtain informa- 

Q_i| tion about the structure of the sources of hadrons in multiple particle 

q" production processes. In this paper the principles of this approach are 

<D . described and some of the controversies about their implementation 

are discussed. 



1 Introduction 

Bose-Einstein correlations among the momenta of identical particles 
produced in a high energy multiple particle production process yield 
information about the structure of the source of hadrons in the pro- 
cess. This has been pointed out in the very first paper on these corre- 
lations Q . In a recent Physics Report issue || Wiedemann and Heinz 
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write: Two-particle correlations provide the only known way to obtain 
directly information about the space-time structure of the source from 
the measured particle momenta. Thus the problem of extracting as 
well as possible the information about the source from the measured 
correlations is of great importance. It is, however, not an easy prob- 
lem. One of the founders of this field of research, G. Goldhaber wrote 
in 1990 Q: What is clear is that we have been working on this effect 
for thirty years. What is not as clear is that we have come much closer 
to a precise understanding of the effect. It is unlikely that the progress 
made during the last ten years would make G. Goldhaber give a much 
more optimistic view. 

There is a great variety of multiple particle production processes. 
In order to illustrate this point we will describe two well-known cases. 
In an electron-positon annihilation at high energy, at first usually a 
quark-antiquark pair is produced. These partons radiate gluons. The 
gluons radiate further gluons, or go over into quark-antiquark pairs. 
After some steps of this cascade, in a process known as hadroniza- 
tion and not well understood, the partons combine into colour-neutral 
hadrons. The hadrons are collimated into two narrow jets pointing 
in opposite directions along the same straight line. A splitting of the 
jets into more jets is also possible. There is an alternative way of 
looking at this process. The first generation quark and antiquark are 
the ends of a colour string. As they fly away the string stretches. 
After some time the string breaks. At the breaking point a quark 
antiquark pair is formed, so that each of the pieces of the string is 
again a string with a quark at one end and an antiquark at the other. 
The string pieces break again and finally short strings appear, which 
go over into hadrons. It is natural to expect that there is a time scale 
r for the hadronization process. Since, however, the system is highly 
relativistic, it is necessary to specify in which frame this time should 
be measured. We choose the centre-of-mass frame and assume that r 
is the longitudinal proper time at the creation of the hadron defined 
by 



T = s/t 2 - Z 2 , (1) 

where t is the centre-or-mass time, when the hadron was created, and 
z is the corresponding coordinate measured along the jet direction (the 
transverse dimensions are less important). An estimate of the velocity 
of the hadron is z/t. These formulae have an interesting implication. 



Hadron production begins at t = r and first slow particles close to the 
interaction point are produced (\z\ small). Only later and further from 
the interaction point (both t and \z\ large) do the fast hadrons appear. 
This production mechanism is known as the inside-outside-cascade @ . 

A very different multiple particle production process are the central 
heavy- nucleus - heavy- nucleus collisions, known also as central heavy 
ion collisions. Here the usual picture is that of two spheres Lorentz- 
contracted into coaxial discs — we consider the centre-of-mass system 
- penetrating through each other. When the discs fly apart, many 
strings are simultaneously stretched in a tube with a transverse radius 
of the order of the radii of the colliding nuclei. For heavy nuclei and 
high energies the strings are so numerous that they merge, e.g. into a 
quark gluon plasma. Then another poorly understood process, known 
as freeze out, converts the plasma (or whatever is the intermediate 
state) into hadrons. 

There are many obvious questions to ask. What is the transverse 
radius of the tubular (?) region, where the hadrons are created? One 
would expect about one fermi or less for e + e~ annihilations and sev- 
eral fermi for heavy ion collisions. What is the formation time t, 
which elapses between the moment of collision and the moment, when 
the last hadron is produced directly? What is the time At between 
the direct production of the first hadron and of the last? There are 
also many model dependent questions. In thermodynamical models 
one asks about the temperature, in hydrodynamic models about the 
velocity of the collective flow etc. 

Let us review now the main results of the famous GGLP paper 
[|]]. Even if one does not quite share G. Goldhaber's opinion quoted 
above, this is certainly a very important paper and some familiarity 
with its content is necessary for any discussion of the Bose-Einstein 
correlations in multiparticle production processes. 

2 The GGLP contribution 

Consider two 7r + -s with momenta pi and P2 produced at ri and T2- 
If the pions were distinguishable, the probability amplitude to ob- 
serve them them both at r would be a product of the single particle 
contributions: 



A D is(r) = exp[z0i + ipi • (ri - r)] * exp[^ 2 + «P2 ■ (r2 - r)] (2) 

In each of the square brackets, <pi is the phase acquired by the particle 
at birth and the other term is the phase accumulated while propagat- 
ing from rj to r with momentum pj. Since, however, the two pions are 
identical bosons, it is mandatory to symmetrize the amplitude and a 
more realistic formula is 



Aund(r) = —7= exp[i(>i + 4> 2 ) + i(pi - p 2 ) ■ r] 

[exp[i(pi • ri + p 2 • r 2 )] + exp[i(pi • r 2 + p 2 • ri)]] . (3) 

Physically this means that we have added coherently the contri- 
bution corresponding to the possibility that the two pions have ex- 
changed their birth points. Classically for different momenta and 
given birth points, if before the exchange the two pions can reach 
the point r, then in general after the exchange they must miss it. One 
should keep in mind, however, that the distance between the points 
ri and r 2 is of the order of a fermi, while the distance between either 
of them and the point r is of the order of a meter. Classically this 
would not help, but quantum-mechanically in this situation the prob- 
ability of reaching r by both pions is in the two cases the same for all 
practical purposes. For comparison with experiment the result should 
be averaged over all the possible pairs of points i"i,r 2 . By averaging 
the amplitude (^) nothing interesting is obtained. Therefore, GGLP 
assumed that one should average the square of the absolute value of 
the amplitude @. Physically this means that the contributions from 
all the pairs of points ri, r 2 add incoherently. Then the probability of 
finding a pair of pions with momenta pi , p 2 is 

< \A Und \ 2 >= 1+ < cos[q • (n - r 2 )] >, (4) 

where q = pi — p 2 and the Dirac brackets < ... > denote averaging, 
with a suitable weight, over all the pairs of points ri, r 2 . Using as an 
example a Gaussian weight function 

p(n,r 2 ) = (2vri? 2 )- 3 exp[-(r? + r 2 2 )/(2R 2 )}, (5) 

where R is a constant with the dimension of length, GGLP found 



<\A Und \ 2 >=l + exp[-q L 2 R 2 }. (6) 

The parameter R can be interpreted as the radius of the sphere, 
where the pions are produced. Thus, finding R from a fit to the 
experimental data yields the size of the hadronization region. Note 
that formulae qualitatively similar to (g) hold for a broad class of 
weight functions. For q 2 = the cosine being averaged equals one 
whatever are the points ri and r%- Thus the right hand side of (0) 
for q 2 = must be equal two. For large values of q 2 , the cosine 
is a rapidly oscillating function of the difference ri — T2- For smooth 
weight functions, therefore, its average is very small and the right hand 
side of (||) equals approximately one. If the weight function contains 
only one parameter with the dimension of length, let us denote it R, 
the width of the region in q 2 over which the right hand side of (^) 
drops from the value 2 to, say, 1.5 must be proportional to R~ 2 for 
purely dimensional reasons. Thus it is not difficult to improve over 
the Gaussian Ansatz. There is a problem, however. 

It is certainly not true that the probability of producing two iden- 
tical pions with momenta pi and P2 depends only on the momentum 
difference q. GGLP could have tried to save the theory by introducing 
some more complicated weight functions, which would depend also on 
the momenta. They found, however, a much simpler and more bril- 
liant solution. Note that in the model the squared modulus of the 
amplitude for distinguishable particles is a constant. Thus one may 
claim that the right hand side of (|6|) is not the two-particle probability 
distribution, but the ratio 

n t \ P(Pi.P2) /->, 

C 2 (pi,p 2 ) = -, r. (7) 

PDMPl,P2j 

Here p(pi,P2) is the experimentally observed two-particle distribu- 
tion and p_Dis(pi,P2) is the distribution, which would be observed, if 
the the identical pions were distinguishable, i.e. if the Bose-Einstein 
symmetrization were switched off. Of course the numerator of this 
expression cannot be measured, which has led to much discussion as 
described in the following section. This is the famous normalization 
problem. GGLP chose again a very simple solution. They substituted 
for pdis the distribution of tt + , tt~ pairs, where of course symmetriza- 
tion does not occur. The experimental results for the ratio Ci could 
be fitted reasonably well with formula (o) and thus both the apparent 



attraction in momentum space among pions of the same charge got 
explained and an estimate of the radius R was given. 

Let us mention one more fruitful idea from the GGLP paper. The 
right-hand side of formula (y) depends on the Lorentz frame, where 
the momenta are measured. Choosing the rest frame of the pair, where 
the energies of the two pions are equal, we can rewrite the result in a 
covariant form: 

C 2 (pi, P2 ) = l + exp[Q 2 i? 2 ], (8) 

where Q 2 is the square of the four- vector p\ — P2- 

3 Normalization 

The problem of finding the best denominator for the function C2 
has attracted much attention. Theorists often suggest to replace 
PDis (Pi; P2) by the product of single particle distributions Pi(pi)pi(p2) 
With this choice, the function C2 becomes the familiar, standard two- 
particle correlation function. Moreover, in many models terms cancel 
between the numerator and the denominator making the formula sim- 
pler. One may also notice that a two-particle density is normalized to 

< n(n — 1) >, while the single particle distribution is normalized to 

< n >. Thus a better choice of the denominator might be to multiply 
the product of single particle distributions by < n(n — 1) > / < n > 2 . 
The advantages and disadvantages of using this factor have been re- 
cently discussed in ||]. Since it does not depend on q 2 , usually it does 
not have much effect on the parameters of the source found from fits. 
The identification of C2 with the standard correlation function, in 
spite of its advantages, is not very popular with experimentalists. In 
order to explain the reason let us consider the following simple exam- 
ple. Consider a high-energy reaction, where the two initial particles 
go over into two well collimated jets. By momentum conservation the 
two jets must be back to back in their centre-of-mass system. Let us 
assume that the orientation of their common axis can point with equal 
probability in any direction — is isotropic. Then the single particle 
distribution of momentum is also isotropic. On the other hand the 
opening angle between the momenta of two particles is either small, 
when the two particles are extracted from the same jet, or close to tt if 
the two particles are from different jets. Thus the correlation function 



exhibits a very large peak for 9 ~ 0. This peak has, of course, nothing 
to do with Bose-Einstein correlations. Experimentalists would pre- 
fer a definition of C2, where the forward peak reflects Bose-Einstein 
correlations and nothing else, because this makes the interpretation 
much easier. Sometimes a compromise is chosen. For instance in e + e~ 
annihilations, where the situation is similar to that from our example, 
the z-axis is often chosen along the jet axis and not along a direction 
fixed in the laboratory frame. With this choice the single particle 
distribution becomes also strongly peaked for small angles and in 
the correlation function the bumps due to the two-jet structure of the 
events are largely eliminated. 

The most popular choices, however, are improvements over the 
choice of GGLP. For instance one uses "mixed" samples, where puis 
is the distribution of pairs of 7r + -ses, but with each tt + taken from a 
different event. Sometimes Monte Carlo generated samples are used, 
with Monte Carlo generators which do not include Bose-Einstein cor- 
relations. This procedure is not very safe, because such generators 
contain a number of free parameters, which are fitted to the data, 
where the Bose-Einstein correlations are present. One also uses ratios 
of functions C<i obtained from the data to functions C2 obtained ac- 
cording to the same prescription from Monte Carlo. The wide variety 
of methods of calculating the denominator of the function C2 is one of 
the reasons, why the comparison of results from different experimental 
groups is very difficult. This is, however only part of the story. Some 
groups correct for final state interactions (mostly coulombic) and/or 
resonances, others do not. Various cuts defining the data samples are 
used. Some groups assume that every negative particle is a ir~ , while 
other have particle identification. Because of all that great care is 
necessary when interpreting the experimental results and their stated 
errors. This difficulty has been known for a long time cf. e.g. H. For 
a more recent (pessimistic) review cf. @. 

4 Beyond spherical symmetry 

The GGLP weight factor is a function of r 2 = x 2 + y 2 + z 2 and, there- 
fore, it is spherically symmetric. It is natural to replace r 2 by an ar- 
bitrary quadratic form in x,y,z, provided the eigenvalues are positive 
so that the weight function can be normalized to unity. Performing 



the averaging of the cosine one obtains 

3 

C 2 (q,K) = 1 + A(K)exp[- ]T R 2 3 (K)q iqj }. (9) 

Here besides abandoning spherical symmetry two improvements 
have been introduced. In agreement with experimental observation 
(cf e.g. ||) a factor < A < 1 has been included and a dependence 
of all the coefficients on K = (pi + p 2 )/2 has been allowed. The 
coefficients denoted R 2 , do not have to be all positive. Out of the 
nine coefficients R 2 j three are eliminated by the symmetry condition 
Rfj = R 2 t . Moreover, choosing the y axis so that K y = and assuming 
reflection symmetry with respect to the x, z plane we have R 2 = 
Ryx = 0- Thus there are four independent coefficients left. The time 
component go can be easily obtained from the identity K^q^ = (p 2 — 
pl)/2 = 0. 

It is convenient to choose the z-axis in the longitudinal direction 
i.e. for central heavy ion collisions along the beam axis and for e + e~ 
annihilations along the jet axis. The y-axis is perpendicular to the 
z-axis and to the vector K. This fixes also the direction of the x-axis. 
With this choice, the x, y, z directions are often referred to as the out 
direction, the side direction and the longitudinal direction respectively. 
The R 2 parameters are denoted R 2 ,Rq, R 2 and R^ || [|. There were 
speculations that the study of R ou t/R s could give clues as to whether 
there is quark-gluon plasma and/or collective flow in the system |y, 
but the results have not been conclusive and a complete study of all 
the parameters seems now to be the best strategy. Other choices of 
parameters are also possible. For instance one can put 

C 2 (q, K) = 1 + A exp[-i^ _ tffi _ R\ q \ _ T 2 q 2 } (10) 
lPH,or 



C 2 (q, K) = 1 + Xexp[-R 2 T q 2 T - R 2 (q 2 L - q 2 ) - (R 2 + R^ 2 ^ - vq 2 )], 

(11) 
where 



VT 



V 



2 



and u = 7(1, Ot)"u). This parametrization proposed in [|ll]] and im- 
proved in |1J] is particularly popular and is often referred to as the 
YKP parametrization. 



5 Time dependence 



In the GGLP picture the production of all the hadrons was instanta- 
neous at some time to- Formally this assumption is difficult to dis- 
prove. Choosing the time to after all the hadrons have been produced 
and interacted and before the time when they were observed, one can 
calculate the distributions at the observation time using the state at 
time to as the initial condition. Whether the hadrons existed before 
time to, is irrelevant for this calculation. In particular, one may as- 
sume that all the hadrons were created at time to- Guessing the initial 
condition at time to is, however, very hard in this approach. There- 
fore, it is more practical to choose a more realistic conjecture about 
the origin of the hadrons, because then the initial conditions are more 
natural and easier to guess. 



In particular, several authors (cf. e.g. |13|, [14|, [15|]) assumed that 
the production of hadrons is from sources, which fly away from each 
other and such that in the rest frame of a source the production of 
hadrons is isotropic. Consider the simplest case of just two sources. 
If their relative velocity is large and the momenta of the hadrons in 
the rest frames of the corresponding sources are moderate, then it 
is very unlikely that two identical pions from different sources have 
momenta close to each other. On the other hand, all the informa- 
tion about the structure of the source comes from pairs with small 
momentum differences |q|. For large values of |q| function C2 is flat 
and carries no information. Consequently, the observed function C2 
contains only information about the single sources and no informa- 
tion about the distance between them. One finds that the effective 
production region is spherical in its rest frame, while the actual pro- 
duction region composed of the two sources is elongated. This is an 
important piece of information. What we observe is not the total size 
of the hadronization region, but the average size of the so-called re- 



gions of homogeneity [16], i.e. regions, where hadrons with similar 
momenta are produced. For this reason the effective hadronization 
regions observed in experiment are approximately spherical, while we 



expect that at high collision energy the actual production regions are 
strongly elongated, because of their stringy origin. 

There is one more interesting result connected with the problem of 
the time span of the hadronization process || . Consider instantaneous 
hadronization from a spherical shell of thickness 5R. The length 5R 
should be reflected in the momentum correlations along the direction 
%out, because the particles, which come from the inner part of the 
shell, are out of phase with those, which come from the outer part of 
the shell. The same effect can be obtained, however, if some particles 
are produced later than others. From the experimental fact that R ou t 
is not particularly large, one concludes that the hadronization process 
does not last very long. This excludes models, where hadrons are 
produced from the quark-gluon plasma in a first order phase transition 
with a large latent heat. This process would be too slow to be made 
consistent with the observed moderate time interval of the hadron 
production. 

6 More quantum physics 

The formulation of the GGLP model is quasiclassical - one talks about 
a pion with momentum p created at point r, which is not quantum 
mechanics. Much work has been done on formulations consistent with 
quantum mechanics. When dealing with incoherent superpositions 
of states one should use density matrices or density operators. The 
starting point may be the density matrix in coordinate representa- 
tion /j(x,x',t) or in momentum representation p(p,p',t). Often it is 
convenient to replace the vectors a, a' by their linear combinations 

a + = -(a + a'); a_ = (a - a') (13) 

In particular, the GGLP results can be derived from the density matrix 

p(p, p', t) = J dr(p|r)p(r)(r|p'}. (14) 

This incidentally shows that in spite of its pseudoclassical formula- 
tion the GGLP model can be translated into respectable quantum 
mechanics. 

The density matrices, however, do not combine explicitly the in- 
formation about the space distribution of sources and the momentum 
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distribution of the final pions. Therefore, in order to derive the space 
distribution of sources from the momentum distributions of the ob- 
served final particles other approaches have been proposed. 

One can use the Wigner function W defined in term of the density 
matrix by the formula 



W(p, 



,x. 



-) 



dp- 



3 *P--x+ 



p(p,p',t). 



(15) 



The properties of the Wigner functions are described in detail in the 



famous review article [17]. For a recent application to the descrip- 
tion of multiple production of identical particles cf. [O]. In a well 



defined sense [17] Wigner's function is the best quantum mechanical 
analogue of the classical phase space distribution. In this formulation 
Heisenb erg's uncertainty principle is easily implemented. The density 
in phase space should not be too large. Quantitatively the condition 

is 



|^( P+ ,p_,i)| 2 <(7rft) 



-3n 



(16) 



where 3n is the dimension of the vectors p± , or equivalently n is the 
number of particles described by the Wigner function W . 

An alternative approach is to introduce the classical position and 
momentum vectors £, it besides the quantum mechanical position and 
momentum vectors x, p. Each pair of vectors f ,7r defines a quantum 
mechanical wave packet e.g. 



(x|e,vr) 




cxp 



1 



X — £) + ITX 



or equivalently in momentum representation 



(pIM 



- exp 

IT 



1 , 



7T 



+ if • (tt - p) 



(17) 



(18) 



Coherent, or incoherent superpositions of such wave packets are legal 
quantum mechanical states. On the other hand, the values of x and p 
are constrained to be close to the values of f and tt. If the distribution 
of the parameters f and tt is calculated from some classical model, the 
result is directly translated into a distribution of the observables x and 
p, which is consistent with quantum mechanics. This approach was 
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pioneered in ref. JLJJ] . For a recent application and detailed discussion 
cf. @. 

Still another approach is to introduce a source function S related 
to the density matrix by the formula 

p(j>,p') = / d A xexp[ip-x+]S(x + ,p + ). (19) 

In this formula x,x',p,p' are fourvectors. The similarity of this for- 
mula to the formula relating Wigner's function to the density matrix 
caused that the source function is often referred to as a Wigner func- 
tion, a kind of Wigner function, a pseudo Wigner function etc. (cf. 
e.g. |22|| , [PUp . In fact the relation between the source function and 
Wigner's function is not unique and may be quite complicated. We 
will not discuss this problem here, but let us note that the density 
matrices or Wigner's function were calculated at a given time, while 
here the time dependence must be known. In simple cases the time 
dependence of a Wigner function can be found and used (cf. e.g. |23[| ), 
but it is quite complicated and unlikely to be found by a simple guess. 
What is more, the source function depends on two time arguments. 

There is one case, however, when the source function has an easy 
interpretation. If the pions are produced by classical currents J, then 



S(x + , P+ ) = J ^^e-^'{J*(x)J(x')), (20) 

where the averaging is over all the incoherent components of the cur- 
rents J(x). 



7 Momentum-position correlations 

In the GGLP model and in many others the momentum spectrum of 
the produced particles does not depend on the position of the source. 
When momentum-position correlations are introduced in the model, 
unexpected results may be obtained. Let us consider for example 
quasiclassical sources with the single particle position-momentum dis- 
tribution 



/j(x, p) = V2-7T-R2 5(p - Ar) exp 



12 



r 2 
2R 2 



(21) 



Correlations of the type p ~ r occur in classical versions of various 
models, e.g. in string models and in models, where hadronization is 
preceded by a rapid flow of hot matter. The average of the cosine 
from the interference term in the present generalized GGLP model is 



(cos[q • (ri - r 2 )]) = A 6 exp 



K 2 " 


exp 


q 2 

4A 2 i? 2 


cos 


"q 2 " 

A 


X 2 R 2 



(22) 



In comparison with the GGLP result for the Gaussian weight function 
p, two changes are striking. A dependence on the sum of momenta 
K appeared and the dependence on the momentum difference q is no 
more Gaussian. For non-Gaussian distributions the effective radius of 
the distribution of sources in space is usually defined by 



R, 



eff 



-(cos[. 



(23) 



d 

,dq 2 / q- 

For the GGLP model with a Gaussian weight function this reproduces 
the usual result R e ff = R- Ln the present case, however, we obtain 

1 



R, 



eff 



2XR 



(24) 



Due to the position-momentum correlations, R e ff becomes inversely 
proportional to R\ 

A more realistic model with similar correlations can be built as 
follows [25], [26]. For the source function we choose 



S(x,p) = 5(t 2 — z 2 — Tq) exp 



(p — Xx) + (p — Ax)_ (p — Ax)y 



xs 



25 2 



24 



2R 2 
(25) 



This formula has a simple physical interpretation. The subscript T 
denotes the vector component transverse with respect to the beam 
axis z. The last term in the exponent implies that the particles are 
created not too far from the beam axis, typical distances being of the 
order of the parameter R. The first two terms in the exponent impose 
the condition p ~ Ax. The subscripts ± refer to the fourvectror 
components ao ± a z . The dispersion of the transverse components is 
of the order of 5j< and that of the longitudinal component of the order 
of <5 2 . The Dirac delta implies that for each drop of the hot matter 
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hadronization takes place after the same longitudinal invariant time 
To- It is useful that for this source function the corresponding density 
matrix in the momentum representation can be calculated in closed 
form. One finds in particular 



|p(K,q) 



6? 



R 2 M 2 \ 
T -i t~^ exp 



8l + R 2 M* T - 2 } 



where Kq is the modified Bessel function, 



exp 



2M| 



l^o(«)| 2 , 
(26) 



a 



M| 



-7q(K t • q T ) 



2 2 2 

4M| v 



-2/2) + 



2iToMrK/rqT 



(27) 
and m 2 ^ = m 2 + p? T and Mj. = X 2 + K 2 ^ are transverse masses of 
the two particles and of the pair. The single particle distribution can 
be calculated from the formula pi(p) = p(p, 0). 

This model was found to reproduce reasonably well the data for 
single particle distributions and for Bose-Einstein correlations in e + e~ 



annihilations at LEP energies 27]. 



8 Final state interactions and (partial) 
coherence 



In the GGLP model and in many later models, after hadronization the 
pions propagate as free particles. In fact we know, that the majority 
of pions is generated in resonance decays. Moreover, there are the 
nonresonant strong and the electromagnetic interactions among the 
pions. All that has been studied for years, but many problems are 
still controversial. 

Strong nonresonant interactions are usually neglected, but it has 



been pointed out [22] that absorption of the produced pions can lead 
to a decrease of the effective radius R e ff with increasing momentum 
of the pair |K|. The effect of resonances seems to be much more 
important. The long lived resonances, mostly r] and rj', usually decay 
far from the centre of the interaction region. Therefore they simulate 
a large hadronization region and produce a narrow peak in the plot of 
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Ci versus q 2 . This peak is narrower than the experimental resolution 
and consequently its main effect is to reduce the parameter A. The 
resonances with life times neither very long, nor very short, like the 
u! resonance, produce for small q 2 a steep rise of C2 with decreasing 
q 2 . Short lived resonances, like the p meson, increase only slightly 
the measured radius of the interaction region. It has been suggested 
pS[|, [J7]] that in e + e~ annihilation it may be difficult to explain why 
after correcting for the resonance effects the Bose-Einstein correlations 
remain as strong as observed in experiment. 

Coulomb interactions can be easily included by replacing in the de- 
scription of the propagation of the pions the plane waves by Coulomb 
wave functions, which leads to the introduction of the so called Gamow 



factor [19]. This is now known to grossly overestimate the effect (for 
a review cf . ]29| ) . The physical reason is that the introduction of the 
Coulomb wave functions describes the evolution of an isolated pair of 
pions, while in reality there are many other pions around, which partly 
screen the Coulomb interaction within the pair. A direct experimen- 
tal argument is that the reasoning leading to the Gamow factor, when 
applied to ■k + tt~ pairs, gives a strong attraction, which is not seen in 
the data. A simple way of correcting for the screening effect is given 
in ref p0[ ] (see also |29| ). A screening radius tq is introduced. The 
interaction potential of the pair is continuous at r = r$ , Coulomb for 
r > tq and constant for r < r$. The parameter ro is chosen so as to 
reproduce correctly the experimental data for it +/ k~ pairs. Coulomb 
corrections calculated in this way are rather small. 

Another assumption of GGLP, which has been put into doubt, 
was that the production process is completely incoherent. Complete 
coherence would kill the effect, but some degree of coherence seems 
difficult to avoid in realistic models. It is easy to write down general 
formulae including partial coherence (cf. e.g. |19[]), but it is not 
clear how to use them fruitfully. At a time it was suggested that the 
parameter A measures the degree of coherence (A = 1 no coherence, 
A = complete coherence), but now it is clear that this parameter 
is strongly affected by the dynamics of the production process, in 
particular by resonance production, and by experimental conditions 
(e.g. particle misidentification) . 
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9 Two recent models 

In order to illustrate how Bose-Einstein correlations are nowadays 
analysed, we present two recent models. The first (cf. M and ref- 
erences quoted there) is based on analogies with hydrodynamics and 
thermodynamics. It is being used to describe central heavy ion colli- 
sions. The starting point is the (single particle) source function 



S(x,p) ~ tut cosh(y — 7]) exp 



rriT cosh y cosh r]T{rr) — Ptj sinh t]t(^t) 



exp 



r\ rj 2 (r - r ) 2 



2R 2 2(Ar/) 2 2(Ar) 2 



In this formula ~ means that a normalization constant has been omit- 
ted, y is the rapidity along the z axis i.e. along the beam direction. 
The pseudorapidity 

1 , t + z , . 

, = -m— . (29) 



The distance from the z axis r^ = \/x 2 + y 2 . The rapidity of the 
transverse flow has been assumed in the form 

rfr(rr) = r lf-^, (30) 

where rjf is a constant. The model contains six free parameters: 
R,T,r]f,Ar],TQ and At. These parameters have been estimated by 
comparison with the data from the NA49 experiment for collisions of 
lead nuclei at an energy of 158 GeV per nucleon. Some of the results 
have interesting physical implications. The parameter R, interpreted 
as the transverse radius of the hadronization region, is about 7 fm. 
From the known size of the lead nucleus one could have expected a 
number about twice smaller. This means that there is substantial 
transverse spreading of the interaction region before hadronization 
takes place. The parameter T interpreted as the local temperature is 
about 130 MeV. This is less than the temperatures found in models 
used to calculate the chemical composition of the produced hadrons, 
which could mean that the spreading is accompanied by cooling. The 
parameter rjf is about 0.35. This is a very reasonable value. The 
velocity of sound in a plasma is about 1/3 (in units of the velocity 
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.(28) 



of light in vacuum). The parameters tq and At are about 9fm and 
1.5fm respectively. It means that the time span of hadronization is 
short compared to the time between the original interaction and the 
onset of hadronization. This would exclude models, where hadroniza- 
tion is a first order phase transition with a large latent heat, because 
in such models the time span of hadronization is large. One should 
keep in mind, however, the authors' warning that the determination 
of the parameter At from the data is poor. This example shows that 
given a model one can extract from the data much interesting informa- 
tion. Little is known, however, about the model dependence of these 
results. 



A very different picture of the Bose-Einstein correlations |3l], 32 
p3j , pHl has been inspired by a string model of the Lund type. This 
model is taylored for e + e~ annihilations. An annihilation is depicted 
as the formation of a string with a quark at one end and an antiquark 
at the other. This string stretches with the speed of light. Then 
somewhere along the string a quark-antiquark pair is produced and the 
string breaks. The pieces stretch and break again. Finally sufficiently 
short bits of strings hadronize. In the z, t plane the trajectories of 
the ends of all these strings form a closed contour. Let us denote the 
area enclosed by this contour by A. The key assumption is that the 
probability amplitude for a given final state is 

Af ~exp[t£A]. (31) 

We have not written explicitly the factor related to the transverse 
momenta of the produced particles. Denoting by 6/2 the imaginary 
part of £ one finds 

\M\ 2 ~ e~ bA . (32) 

This result is well known from the Lund model. For the description of 
the Bose-Einstein correlations, however, it is the real part of £ which 
is the important one. It is expected to be of the order of the string 
tension i.e. of the order of 1 GeV/fm. The point is that the produc- 
tion amplitude has to be symmetrized with respect to exchanges of 
identical particles. Such an exchange, however, changes the area A 
and because of the nonzero real part of £ the phase of the amplitude. 
Numerical calculations show that this model gives a reasonable de- 
scription of the Bose-Einsten correlations in e + e~ annihilations. This 
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is very interesting, because this model, contrary to all the previous 
ones, does not contain incoherent components, random phases etc. It 
is curious what happens, when in a process more than one string is ini- 
tially produced. In e + e~ annihilations there are events, where two W 
bosons are produced, which implies two strings. It has been suggested 
pi] that in this case only pions from a single string should exhibit 
Bose-Einstein correlations. To be sure, nobody doubts that pions are 
bosons and that consequently their production amplitude should be 
suitably symmetrized, but the attraction in momentum space, known 
since the GGLP paper as Bose-Einstein correlations, requires in addi- 
tion certain phase relationships, which may not be realized for pions 
originating from different strings. In central heavy ion collisions many 
strings are produced. By extension of the previous argument one could 
expect very weak Bose-Einstein correlations, which experimentally is 
not the case. This suggests that the state from which hadronization 
occurs in heavy ion collisions is not a bunch of strings. The strings 
must somehow merge and form a very different object, perhaps a vol- 
ume of quark-gluon plasma. 

10 Conclusion 

Bose-Einstein correlations attract much interest. Many hundreds of 
papers have been published on this subject. In principle they offer the 
only access to some important information about the hadronization 
process. Typical questions are: what is the size and shape of the 
hadronization region, how is hadronization distributed in time, what 
is the coherence of the sources, what is the hadronization mechanism 
etc. In practice it is difficult to find model independent, definitive 
answers to these questions. The common approach is to study models, 
which either demonstrates the viability of certain scenarios, like the 
sting model discussed above, or give tentative answers to the questions 
concerning hadronization. 
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